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Does the gamma-ray flux of the blazar 3C 454.3 vary on 
sub-hour time scales? 
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ABSTRACT 

In the early days of April 2010, the blazar 3C 454.3 (z = 0.859) underwent a strong 
7— ray outburst, reaching fluxes {E > 100 MeV) in excess of 10^^ ph cm~^ s~^. 
The Fermi Gamma ray Space Telescope performed a 200 ks long pointed observation 
starting from 5 April 2010 19:38 UTC. This allowed us to try probing the variability 
of the 7— ray emission on time scales of hours or less. We found the variability on 
a few hours time scale. On sub-hour time scale we found no evidence of significant 
variability, although the present statistics is not yet conclusive and further observations 
are needed. 
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1 INTRODUCTION 

Blazars are one type of active galactic nuclei (AGN) with 
relativistic jet. The small viewing angle of the jet makes it 
possible to observe strong effects of the special relativity, 
such as a boosting of the emitted power and a shortening 
of the characteristic time scales. Although the variability 
in blazars has been observed at all the frequencies (see, for 
example, the reviews by Ulrich et al. 1997 and Wagner 2008), 
the part at high energies deserves a specific interest, since 
most of the bolometric power of these sources is emitted at 
7 rays. 

The shortest time scales measured by the EGRET ex- 
periment onboard the Compton Gamma-Ray Observatory 
are ~ 4 hours for PKS 1622-29 (Mattox et al. 1997) and 8 
hours for 3C 279 (Wherle et al. 1998). Early results from the 
Large Area Telescope (LAT) onboard the Fermi Gamma-ray 
Space Telescope (hereafter Fermi) indicated similar results: 
from about half-day for PKS 1454-354 (Abdo et al. 2009a) 
and PKS 1502-1-106 (Abdo et al. 2010a) to 5 - 6 hours in 
the cases of 3C 454.3 and PKS B1510-089 (Tavecchio et al. 
2010) and 3C 273 (Abdo et al. 2010b). 

These time scales are in agreement with the common 
paradigm that the characteristic spatial scale of the emitting 
region is of the order of the gravitational radius rg = GM/c^ 
of the central spacetime singularity (Begelman et al. 2008). 
The zone where most of the dissipation occurs is located at 
distances greater than T^rg, where F is the bulk Lorentz fac- 
tor (for most blazars F ~ 10, see Ghisellini et al. 2010). The 
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size of the emitting region is in turn linked to the variability 
time tvar by the relationship R < cfvar'5(l -I- z) (where 5 is 
the Doppler factor). Obviously, it is expected that R > rg. 
Moreover, the emitting region must be optically thin and 
sufficiently far from the central black hole to allow the 7 
rays to escape, otherwise they convert into electron-positron 
pairs. 

An order-of-magnitude estimation of the above param- 
eters results in an overall agreement with the observed vari- 
ability on scales of days-hours. However, the recent detection 
by ground-based Cerenkov telescopes (HESS and MAGIC) 
of fast variability (minutes time scale) at 7 rays in some 
TeV BL Lac Objects severely threatens the above scenario. 
In 2006, during an exceptional outburst (average flux 7 Crab 
with peaks of more than 14 Crab), PKS 2155-304 displayed 
variability with doubling flux time scale of about 200 s (Aha- 
ronian et al. 2007). In 2005, Mkn 501 changed its flux within 
a few minutes (Albert et al. 2007). Given the masses of these 
two blazars of the order of W^A4q, the measured variability 
is more than one order of magnitude smaller than the mini- 
mum allowed. Several solutions have been proposed, from a 
"simple" increase of the Doppler factor {5 > 50) to invok- 
ing an internal structure of the jet (Begelman et al. 2008, 
Ghisellini & Tavecchio 2008, Giannios et al. 2009). 

As noted by Begelman et al. (2008), while these expla- 
nations could flt reasonably with BL Lac subclass of blazars, 
the same is not appropriate for flat-spectrum radio quasars 
(FSRQ) subclass. Indeed, this type of blazars generate 7 
rays by inverse Compton on a population of seed photons 
external to the jet (external Compton, EC). Therefore, the 
pair photosphere can be at fairly large distances (> 10*rg), 
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which in turn means that the blob has a size so large to de- 
termine an insufficient energy density to develop strong and 
fast flares. 

It is evident that the search for sub-hour time scales in 
FSRQ can have strong impact on the current knowledge. If 
measured, such very short variability would call for a strong 
revision of the present day models and understanding of 
relativistic jets. What we need is a strong outburst - like 
those occurred in the BL Lac Objects PKS 2155-304 and 
Mkn 501 - and a highly performing instrument. 

To date, the shortest time scale observed at high ener- 
gies in a FSRQ is ~ 2000 s, but in the hard X rays. It is the 
case of NRAO 530 as observed by INTEGRAL (20-40 keV 
energy band) during an exceptional flare in February 2004 
(Foschini et al. 2006). Despite of the good performance of 
INTEGRAL, that episode consists of a single 5(t detection, 
while the source is below the detection limit of the instru- 
ment during the remaining of the time. No multiwavelength 
observations were possible at that time and, speciflcally, no 
7— ray satellites were operating. Therefore, this remains an 
isolated exceptional case and, as for every exceptional claim 
without a strong observational support, there is always an 
"halo" of doubts. 

The launch of Fermi /LAT in June 2008 could give more 
opportunities to search for sub-hour variability. The Large 
Area Telescope (LAT, Atwood et al. 2009) could have the 
necessary sensitivity to probe the shortest time scales at 7 
rays in quasars. If we consider a flux of 10~^ ph cm~^ s~^ 
{E > 100 MeV) and an eflective area of ~ 5000 cm^ (see 
Fig. 14 of Atwood et al. 2009, by taking into account the 
0.1 — 100 GeV energy range), it results that for normal in- 
cidence there are about 15 photons per 5 minutes (about 
a 4(T detection, by assuming that on such short time scales 
the background is almost absent). However, LAT operates 
almost always in scanning mode. This, in turn, has the great 
advantage to offer an efficient monitoring of all the sky (a 
complete coverage every three hours or two orbits), but it 
hampers the possibility to probe time scales shorter than 
three hours (the source is not always at the boresight of the 
instrument and, therefore, the effective area is smaller). 

A first tentative to bypass this issue has been recently 
made by performing a pointed observation during the giant 
outburst of 3C 454.3 (z = 0.859) occurred in early April 
2010 (Wallace et al. 2010), when the source reached fluxes in 
excess of 10"^ ph cm"^ s"^ {E > 100 MeV). Such flux levels 
were already reached by 3C 454.3 during the first ten days 
of December 2009 (see Bonnoli et al. 2010 and references 
therein), but it was not possible at that time to perform a 
pointed observation and hence to search for sub-hours time 
scales. 

Here we report the study of the variability of the 7— ray 
emission of 3C 454.3 as observed by Fermi /LAT during the 
pointed observation performed from 5 to 8 April 2010. 

2 DATA ANALYSIS 

Data of Fermi /LAT (Atwood et al. 2009) have been down- 
loaded from the publicly accessible web site of HEASARCQ 
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and analyzed by means of the public software package LAT 
Science Tools v. 9.15.2, including the calibration files 
(Instrument Response File, IRF P6_V3_DIFFUSE), the maps 
of the spatial diffuse background and the spectrum of the 
isotropic (cosmic and instrumental) background. 

We retrieved the photon list with E > 100 MeV and 
spacecraft: data from 5 Aprfl 2010 19:38 UTC to 8 Aprfl 2010 
03:12 UTC (MJD 55291.8180-55294.1333) corresponding to 
an elapsed time of about 200 ks (see Fig. [T} . Events of class 
3 (diffuse) inside a radius of 7° centered on the coordinates 
of 3C 454.3 {RA = 22'^ : 53"" : 57".7; Dec = +16° : 08' : 53", 
J2000) and coming within 105° of the zenith angle have been 
selected. This data set has been further filtered to exclude 
the events outside the good-time intervals and within 8° 
of the Earth limb. Livetime and exposure maps have been 
calculated, as well as the diffuse background. 

An unbinned likelihood algorithm (Cash 1979, Mattox 
et al. 1996), implemented in the gtlike task, has been used 
to extract the fiux and the photon index of the source emis- 
sion. A single power-law model in the form F{E) oc E~^ , 
where F is the photon index, was used to fit the energy fiux 
distribution of the blazar. It is known that the 7— ray spec- 
trum of 3C 454.3 displays a break at a few GeV (Abdo et al. 
2009b, Finke & Dermer 2010), but since we want to probe 
the shortest time scales (hence, lower statistics), it is better 
to adopt a single power-law model. The broken power-law 
model has been used to fit the data of the whole observation 
(see the next Section). The backgrounds (spatial diffuse and 
isotropic) measured during this fit are kept fixed during the 
processing of the light curves with hours time bins or less, 
since we do not expect strong changes of the background in 
such short time scales. 

A first run with the optimizer DRMNFB is done to grossly 
estimate the spectral parameters. Then, the output of this 
process is used as input for a second run with the NEWMINUIT 
optimizer, to better assess the measurement errors. In the 
following, only statistical errors are displayed. The use of 
the P6_V3_DIFFUSE response file guarantees that systematics 
affect the results by a value of less than 20% (Rando et 
al. 2009). However, since this work is based only on LAT 
data, with no comparison with other instruments at other 
wavelengths, there is no need of absolute fiux calibration. 
Hence, systematics are not applied. 

This procedure is performed iteratively for each bin 
time within each good-time interval in order to generate the 
light curve. We discarded all the detections that have no 
sufficiently high Test Statistic {TS, see Mattox et al. 1996 
for the definition), i.e. we considered only time bins with 
detection at a ~ VtS > 5. A light curve with 3 hours time 
bin is displayed in Fig. [U left panel. The effective source 
ontime changes within the elapsed time bin, ranging from 
10 to 98% in the 30-minutes bin light curve. Therefore, we 
checked for a reasonable high TS even in presence of low 
source ontime. For example, the bin with the lowest source 
ontime (10% of the elapsed time bin) has TS = 73. 

We note that the average TS increased by a factor ~ 
2.5 during the pointed observation, as expected from the 
increase of the instrument sensitivity when changing from 
scanning to pointed mode (the light curve in scanning mode 
is not shown). 
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Figure 1. Light curves of 3C 454.3 {E > 100 MeV) with 3 hours (left panel) and 1 hour (right panel), respectively. Time starts on 31 
March 2010 (MJD 55286), so that the number of the days corresponds also to days of April. The horizontal dashed lines in each panel 
correspond to the weighted average plus/minus one standard deviation. 



3 DISCUSSION 

The whole data set of 200 ks observation has been fitted 
with a broken power-law model in the form of F{E) oc 
(S/Sb.eak)''^ for E < Streak and F(E) oc (-E/^break)''" Oth- 
erwise. This results in these parameters: the photon indices 
below and above the break energy are Fi = 2.22 ± 0.04 
and r2 = 2.80 ± 0.08, respectively. The break energy is 
0.9 ±0.2 GeV and the flux integrated over the 0.1-100 GeV 
energy band is (12.7±0.3) x lO"*^ ph cm s . The recorded 
counts are 3044 and TS = 10592 (> lOOcr). We compare 
these values with the early observation with Fermi referring 
to the first month of data (Abdo et al. 2009b). In August 
2008, the flux was about a factor 4 smaller than in April 
2010, with a value of (3.0 ± 0.1) x lO"'' ph cm"^ s"S the 
break energy was greater (2.4 ± 0.3 GeV), and the pho- 
ton index above the break quite softer, with a value of 
r2 = 3.5 ± 0.2. Fi is consistent within the measurement 
errors. It seems that the higher flux needs of a harder spec- 
trum. 

The 7— ray emission of 3C 454.3 in Fig. [T] (left panel) is 
variable on a few hours time scale. A fit with a constant-flux 
line gave X^/dof = ~ 2.7 (with 18 de grees of freedom, 
this means a chance probability of 0.99). However, we note 
that the last point, which is mainly driving the global vari- 
ability of the light curve, is not a statistical fluctuation, but 
the beginning of a real flux increase. This is confirmed by 
the fact that on 8 April there was the peak of the 7— ray 
emission of the April 2010 outburst, with a day average of 
~ 16 X lO"'' ph cm s . On the other hand, no strong 
variability of the photon index has been measured, although 
some trend is visible between ~ 6.5 and ~ 7 April: F changed 
from 2.58 ±0.13 on ~ 6.63 April to 2.21 ±0.08 (~ 3cr signif- 
icance) on ~ 6.88 April. Similar spectral changes have been 
observed in 7— ray blazars by EGRET (Nandikotkur et al. 
2007) and LAT in the case of 3C 273 (Abdo et al. 2010b). 

The light curve with 1 hour time bin (Fig.[T] right panel) 
has larger error bars and hence the slight variability seen 
in the 3-hours-bin curve is now smoothed. The fit with a 
constant flux line gives = 1-6, but the visible trends are 
consistent with those observed in Fig. [T] (left panel). The 
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Figure 2. (top panel) Flux light curve of 3C 454.3 {E > 
100 MeV) with 30 minutes time bin; (bottom panel) photon index. 
Time starts on 31 March 2010 (MJD 55286), so that the num- 
ber of the days corresponds also to days of April. The horizontal 
dashed lines in each panel correspond to the weighted average 
plus/minus one standard deviation. 



most significant trend is the drop in flux observed between 
~ 6.7 - 6.8 April. 

To evaluate the time scales of these trends, we calcu- 
lated the time of exponential rise or decay as defined by: 

F(t) = F(to)exp[-(t-to)/r] (1) 

where F{t) and F{to) are the fluxes at the time t and 
to, respectively, and r is the characteristic time scale. 

In this case, the drop on ~ 6.7—6.8 April visible in Fig.[T] 
(right panel), has r = 4.8±4.6 hours (significance of the flux 
variation 3.1a). This value is confirmed in the 30-minutes 
time bin light curve shown in Fig. [2] (r = 4.7 ± 4.4 hrs, 
S.lcr). Moreover, the rise at the end of the curve between 
~ 8.0 - 8.1 has r = 2.7 ± 1.0 hrs {3.3a). 

It is possible to estimate the minimum Doppler factor S 
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(e.g. Dondi & Ghisellini 1995, Mattox et al. 19970, although 
during this pointed observation it was not possible to per- 
form multiwavelength observations, because the blazar was 
too close to the Sun. Therefore, we adopted for the X-ray 
flux and spectrum, the values measured during the Decem- 
ber 2009 outburst, when 3C 454.3 reached similar 7— ray 
fluxes: Q = 0.4 and -Fi kcv ~ 30 /iJy (Bonnoli et al. 2010). 
We obtain 5 ^ 14, for 1 GeV 7 rays and r = 2.7 hours, a 
high value for a lower limit, but not unlikely (cf Ghisellini 
et aL 2010). 

We noted a few cases of sudden drop in the flux with 
~ 30 min time scale on ~ 6.1 — 6.2 and ~ 7.3 April. A closer 
inspection of the data revealed that these bins had small 
source on time (~ 25% of the whole bin). Therefore, it is 
likely that it is a fake drop, caused by the not good recon- 
struction of the flux in presence of too few events, although 
sufficient for a high TS. 



4 FINAL REMARKS 

We presented the variability analysis of the 7— ray observa- 
tion with Fermi /L AT of the blazar 3C 454.3 pointed dur- 
ing the outburst of early April 2010. Contrary to the usual 
scanning mode, this time Fermi performed a 200 ks pointed 
observation (Wallace et al. 2010), which guarantees a higher 
sensitivity and a longer source ontime. These particular set- 
tings, together with the exceptionally high flux emitted by 
the blazar, enabled us to probe the shortest time scales that 
we found to be in the range of ~ 2 — 5 hours, confirming the 
indications found by Bonnoli et al. (2010), Tavecchio et al. 
(2010) and in agreement with the current theories. 

Although the significance of the events is not very 
strong (~ 3a), the fact that it is present in the light curves 
with different time bins is an indication that the findings 
are well grounded. On the other side, we have not found 
significant evidence of sub-hour time scales, but the present 
data set is not yet conclusive. We noted that the fiux in- 
creased at the end of the pointed observations with r ~ 2.7 
hours. Therefore, we could expect (or reasonably hope) that 
hours or less time scales may be effectively probed with 
Fermi /LAT under conditions of even greater fiuxes, likely 
as occurred in December 2009 (Bonnoli et al. 2010). 
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^ In the present work we used the most recent value for the 
Hubble-Lemaitre constant Hq = 73 km s~^ Mpc~^ (Freedman 
& Madore 2010). 



